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ABSTRACT 


This thesis is an accuracy study of three empirical 
series system lower confidence limit procedures, Computer 
Simulations were used to compare the accuracy of the pro- 
cedures using the same set of data with 500 replications 
for each case. Modifications were then made to improve the 
accuracy of the Log-Gamma method, 

The systems simulated had reliabilities ranging from .6 
to .95. They were composed of four, fifteen, or forty com- 
ponents, and had component sample sizes of ten, twenty, or 
EDU CY, 

The Lieberman-Ross method was accurate but had high 
variance. The Log-Gamma method was conservative with small 
Sample sizes, but became more accurate with increased test- 
ing. Randomization of the Log-Gamma method was the most 


successful of the modifications attempted, 
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I. INTRODUCTION 


There exist many lower confidence limit procedures for 
Series system reliability which are based on component data. 
All of them are approximate in one way or another for real- 
istic testing situations, Two such procedures are the 
Lieberman-Ross (LR) procedure and the Log-Gamma (LG) proce- 
dure. The LR lower confidence limit (ICL) is based on 
operating time data and assumes time to failure is expo- 
nentially distributed. The LG ICL is based on attributes 
data and thus is non-parametric. A modification of the LG 
procedure, called the Continuous Log-Gamma (CLG) procedure, 
is based on time data and assumes an exponential failure 
distribution. 

In this thesis the accuracy of the LR, LG, and CLG pro- 
cedures are compared when data is generated from the expo- 
nential distribution for serial systems with four, fifteen, 
and forty components, Sample sizes are generated ina 
manner which is somewhat representative of data accumulated 
from missile flight tests. 

The accuracy analysis was extended to the case where 
all components were assumed to have a Weibull distribution 
with increasing and decreasing failure rates, 

In all cases, the LR and LG procedures were slightly 
modified from their more common forms in order to accommodate 


zero failures on all components. 





Since the component test data was generated using the 
exponential assumption, the LR method enjoyed an advantage 
over the LG method in that continuous WL procedures are 
exact while discrete LCL procedures are not. All discrete 
100(1-a)% ICL's, Rs (a)? for system reliability R. have 


the property that 
A 
PIRS uta) SR J 21-0 


For some values of R., this probability is 100%. Thus dis- 
crete ICL's will usually be lower in value than a continu- 
ous ICL using the same set of data. Discrete LCL's have 
the advantage, however, of remaining the same regardless 

of the underlying distribution. Consequently, they are 

not susceptible to inaccuracies due to a mistaken assump- 
tion about the probability distribution of the operating 
component time data as are continuous procedures, That is, 


A 
3 Rs L(a) is discrete 


P[R <R,J21-a 


S,L(a) 


regardless of the probability distribution of the component 


A & . 
time data, whereas if Rs L(a) is continuous 
[Re <R | = 1 


provided the required assumption made about the probability 
distribution of the component failure data is correct, 


Even if the assumed continuous distribution is not correct, 


10 





however, i M may Still be reasonably accurate, depend- 
ing upon the true component failure distribution. Of course 
these distributions are never known. 

Finally, continuous IOL's required more refined data 
than discrete ICL's, Discrete LOL procedures can be used 
if continuous time data is available but continuous ICL 


procedures cannot be employed with only attributes data. 


IM 





II. DESCRIPTION OF THE METHODS 


A. THE LOG-GAMMA METHOD 

The following development is a summary of that given 
in Ref, 1. The Log-Gamma (LG) method requires only attri- 
butes data, and assumes nothing about the failure distribu- 
tion. It is designed for independent series systems, does 
not require equal component sample sizes, and is approxi- 
mate. A modification of the LG method for systems composed 
of a mixture of series and parallel-connected components 
was explored by Ref, 2. 

Suppose a system consists of K components in Series, 
and that N; copies of component i are put on test and oper- 
ated until failure or completion of mission, whichever 
occurs first, Denote F. as the number of components of 
type í that did not complete the mission. Then the maximum 


likelihood estimate for system reliability is 


A K A 
R = Tí R. (igi) 
Se Ti 
N, =- F 
where R Ze (2,2) 
i N, 


In the LG procedure the method of moments is used to 
A 
fit the random variable -inR, with the two-parameter gamma 


distribution. The gamma is then transformed into a chi 


12 


7 





ove 
a 
FF ~~ 


MEM s ہمیچ جعجحجہے‎ 


square distribution, about which probability statements 
are made and a lower confidence limit (LCL) obtained. 


That is, define 


K K 
S = -1nR, = -1n TT R= > È in(1-Q, ) (#2 3) 
i=1 = 
K > K (See Appendix A) 
x © (Q,+3Q,)= ET, 
1 = (2.4) 
A 
An unbiased estimator Ti for Ti is 
22 
۲ 9۵ . pt (2.5) 
A 19; u 2 „5 
21 = 3 
where As = EN - 2 (2.6) 
N; ۱ 
B; = MT (2.7) 
A و‎ . : 
and a Se = 1756 (See Appendix B) 
i Ny 
A 
That Ti is unbiased is important, because 
^ K ^ 
S= LIT. (2.8) 
i-i + 


is used as an estimator for S, thereby accumulating any 


^ 
bias present in the Ti. Àn approximate value for the vari- 


^ 
ance of S is 


^ ^ 
Var(S) - Var(T 


po 
I MA 
a 


ið 


a 
۰ مھ‎ 


(2.9) 


22 
pe 
1 P5 
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Appendix C derives the actual variance of A and shows 


A 
equation (2.9) is an excellent approximation for Var(S). 


Next, fit S with a gamma distribution. The probability 


distribution of $ is then given by the density function 


1 EJ 
fa(xir, 0) NT Exp( -x/ 0) x>0,r>0 (2.10) 
4> 0 
It follows that 
E(S) = re (2,11) 
Var(S) = a“ 
A 
Since S is unbiased, 
E(S) =S 
K 
= I T; (2.13) 
31 


Solving equations (2.9), (2.11), (2.12), and (2.13) 


Simultaneously gives the shape parameter 


K Æ Ti " 
` TT L ت‎ 2l 
isi ? z N; 


and the scale parameter 


K 2 K 
i=1 N; 


Thus, r can be estimated by 


. KNEE 
r= ALT, Ë FT m (2,16) 
i-1 i hel Y 
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(2,12) 


e 





A 
== 





Since $ is distributed gamma (r,ð), is distributed E. 
Then 


۳ 
l-a = P< > UT. a, 21? (2. 17) 


2 
P(2r$2 r6 «Ji 4 2.) 
9 


Since S is unbiased 


E(S) = S 


= -lnR (2,18) 


Equation (2.17) becomes 


) 


A 2 
1-a = P(2r8>-10R, ٦ 


=, 2r 


^ 
-2r3 
2 <R,) (2.19) 
1-0,2r 


Therefore, the LG 100(1-a)% ICL for system reliability Ra 


P(Ex 


is 
-2%8 
Ro ETE Exp 1 1 (2.20) 
leg.2r 
To preclude usage of non-integer degrees of freedom, the 
approximation 


„el (2,21) 
dmn [or] 


is used in equation (2.20), where Zr] denotes the smallest 
integer greater than or equal to 2f. Approximation (2.21) 
‘was shown to have little effect on the IG procedure accur- 


acy (Ref. ھا ہے‎ 13]. However, it may produce monotonicity 


15 





problems under certain circumstances (discussed in Chapter 


V). Thus, the 100(1-a)% LCL becomes 


A 


Rs, Ila) ~ rn P 


and was the LG LCL initially simulated by this study (Chap- 
ter IV). 


B. THE CONTINUOUS LOG-GAMMA METHOD 

The Continuous Log-Gamma (CLG) method is developed in 
detail by Ref. 3, pages C-1 through C-20 (as the "Classi- 
cal” method). It differs from the LG method only in that 
it requires component time data and assumes exponential 
failures, 

The derivation of the CLG procedure is analogous to 
that of the LG procedure, and is outlined as follows, 


Define 
Ze -InR, 


一 R. (2,23) 
i=1 Í 


Using time data and assuming exponentiality, Ri can be esti- 
Matec by 


A 
Ri = Exp( = ) iz1,K 


Tsum; 
where Tsum; is the accumulated operating time of the ith 
component., Then Z is estimated by 


A K p 
7 E MMC NE (2,246) 
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The method of moments gamma fit to 2 uses 


^ K F; 
= Ms Tess (2.25) 
A K F; 
Var(Z) = D Tsun, Ni (2.26) 
3(1 
8302 2ص‎ (2,27) 
and Var(Z) = LA“ (2728) 


where L is the shape parameter and Á is the scale parameter, 


The resultant CLG 100(1-a)% ICL is 


Re (0) ees (gli (2.29) 
-q 


h ÁR be Vake | | 
; =1 Tsum, A (res l 


To preclude non-integer degrees of freedom 
2 a: | 
1-0, [21 |) — 


where [at | is the smallest integer greater than or equal 
A 
to 2L. 


C. THE LIEBERMAN-ROSS METHOD 
1, The Original Procedure 
The Lieberman-Ross (LR) method was designed as an 
exact lower confidence limit (LCL) procedure for independ- 
ent series systems (see Ref. 4). It requires component 
times between successive failures from an assumed exponen- 
tial failure distribution, and does not require equai com- 


ponent sample sizes, 
17 








The test plan may be described as follows, Suppose 
N; copies of each component type i are available for testing, 
Place a copy of the ith component on test, and replace it 
when it fails. Continue testing with replacement until all 
available copies of one component are exhausted (this proce- 
dure need not be conducted simultaneously for all component 
types). For the ith component, denote the time to the first 
failure as Tias between the first and second failure as Ti2) 
and so forth, In general, the time between the jth and 


B th failure on the ith conponent is T.,, i-21,.....K 


1) 


క Ni 


where the Ti are in mission units, Reference 4, p. 840, 
also describes the LR procedure for Type II censoring at 
the rth failure of R items on simultaneous test, Reference 
3, p. C-32, illustrates the procedure using component times 
between successive failures and the observed success time 


of the last copy of each type tested. 
The LR ICL is computed as follows, Define 


U ale 0 (2.32) 


Let I; be the greatest index such that 


I; 
E رر‎ 1 నా (2.33) 


Also define 
(2.34) 


H 
u t3 7 
I 


i-i ! 


18 


< 





The variable I is the total number of failures in all con- 
ponents up to and including U (see Figure 1). Then the LR 
100(1-a)% LCL for system reliability R. is given by 


Re MES ) (2555) 


Figure 1 
Depiction of Lieberman-Ross Parameters 
U=35 and I=6, when K=2, N1 =+， and N,=3 


COMPONENT 


i=1,N,=4 


= ——— 4f 


5 10 15 20 25 30 35 40 


U Time 


To implement tne LR procedure, the component data 
order must be specified, If only a set of component data 
were available, then a different procedure would have to be 
used, Note that the opposite chi square tail is used from 
that of the LG and CLG methods, 

2. The Modified Lieberman-Ross Method 





The original procedure utilizes accumulated compo- 
nent times between failures, The Modified Lieberman-Ross 
procedure (MLR), uses instead component times to one mission 
unit or failure, whichever occurs first. The MLR metnod is 


thereby more readily applied than the LR method to actual 


19 


system test data, especially when testing is destructive 
(for example, a missile firing). The MLR method also can 
be used when no failures occur, whereas the LR method could 
not (see Appendix D). 

The MLR procedure is 

a. suppose N; copies of the ith component are avail- 
able for testing. Define Tij as for the LR procedure, ex- 
cept all times past one mission unit are truncated to one, 


Next, define 


* 14 where i4 is the index of 
In. = LOT... the first component time E 
i ju lJ less than one mission (2.36) 


unit for the ith component 


S io where i, is the index of 
Te = Y کت‎ the second component time 
12 ل‎ 2+1 less than one mission unit (2.37) 
for the ith component 


where iF: is the index of 
m" İF; the last component time 

iF; = L T; y less than one mission unit (2.38) 
7 ت‎ for the ith component 
3ئ‎ 


Then obtain 


F: 

+ 1 * 

T; = SS T: ; 191... NS (2.39) 
j=1 


b. Compute the minimum accumulated times between 
successive failures in any component type. 


K 


- { min m” s > 
E | Ti a p 
| nin N; if Ë F;=0 (see Appendix D) 
== i=l 
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4 
The variable I. is now defined in terms of the Ti. 


the greatest index such that 


h = 
1 I; m 
E I اع‎ if F.>0 (2,41) 
: 4 ii) = i 
j=1 
0 1Í F; =0 
Then 
K 
1 2 9۵52 و41‎ I I, (2,42) 
i=1 ` 
c. Now the MLR 100(1-a)% LCL for Ro can be 
computed 


R m E (2,43) 
Sr (air 2) 20 m 


Like the LR method, the component data order must be 


Specified. 
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IN THE SIMULATION 


A. GENERAL METHODOLOGY 
Any lower confidence limit (LCL) procedure can be evalu- 
ated by computer simulation of its distribution as follows: 
و1‎ Suppose it is desired to evaluate a proposed 
100(1-a)% LCL procedure, denoted Eun for system relia- 
(+4 ۸.0 Roe Then the assertion is 


P[Rs 10) £8s 2 ۲ 0 (3.1) 


Equality should hold if Re TOS is a continuous random 
۱ ` 
variable, 
ga Rg = f(R,R,... 4, RJ, where Ri is the true reliabil- 


ity of the ith component, For an independent series system, 


K 
R„ = TT Ri. If exponentiality is assumed, an equivalent 
0 


statement is Re = fO sòns eo... sA) e 
3. Assign values to the parameters a, K, and R, (or Ai), 
i21,2,....K. Simulate life testing by generating random 


component times-to-failure, and compute the resultant 
A 


BENE) j A 
4, Generate the approximate distribution of Ro L(a) 
9 


by repeating step 3 five hundred times. 


A 人 
5. Order the Rs La) realizations to get o” 


... E 0) అ 
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6, Find the (1-q)500th order statistic of Ra po) and 
denote it A, ,. Thus A, , is the (1-a)th percentile of 
the distribution of the LCL random variable Re ra) If 
in fact, Es La) is an exact 100(1-a)% ICL procedure for 
Ro» then A should be approximately equal to R 

7. By repeating steps 3 through 6 for various sets of 
K,NioNoro... No اوه همه وحو و‎ and a, and comparing the 
resultant ۸1 to R the overall performance of the WL 
procedure can be evaluated, 

The actual confidence level given by رت وا‎ can te ob- 
tained by finding the order statistic of the generated 
distribution which matches (or closest matches) Re! 
the index of this order statistic is denoted కూ then 


人 


۰ * 
| 
(100%) 255 = Actual Level of Confidence of Rs (ల) iene) 


Equivalently, if తడ Ego the procedure is a conservative 
one, and vice-versa, 

Figure 2 depicts a conservative outcome with A 1 uq E 
The degree to which Maa equals ie reflects the validity 


of equation (3.1) and therefore the accuracy of Ro ND 
9 


Figure 2 


A 
SewL(a) 
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Not only should an LOL procedure be accurate, but also 


A 
should have small standard deviation. Even though R 


Se) 
may be an exact LCL procedure, its practical value is 
diminished if it produces widely dispersed ICL's, This 


study used both criteria for comparison of the three methods. 


B. GENERATION OF SAMPLE SIZES 

The component test data are generated in a manner some- 
what representative of data accumulated in missile flight 
tests. Suppose some number of Series systems are to be 
used for testing. Failure of a component may fail the sys- 
tem, but still allow testing on components upstream of the 
failed one. For example, if the system is a five-stage 
missile, failure of the third stage may or may not allow 
testing on stages four and five, 

Input parameter N,, the number of mission tests on the 
first component, is the number of systems to be tested, A 
component failure is assumed to preclude testing on upstream 
components one-half the time, The simulation computes the 
number of tests on the ith component as 

i-1 
N; = Ny ar a | (353) 
where అ ۰ denotes "the smailest integer greater than or 
equal to", To illustrate, let N,=20 and assume the first 
component had three failures, Then the component tíns 


Sample size becomes 


N, = 20 -[3(3)| = 18 mission trials 
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It should be understood, however, that the three procedures 
themselves do not require this scheme, but are valid for 


varying sample sizes, 


౮, THE SIMULATION ALGORITHM 

Step ones Initilize the number of series components K, 
the number of systems to be tested N the true reliability 
of the ith component Res Y-cl2,.4,h, and the number of 
replications desired ("Do Loop” index). 

Step two: Using the probability integral transform 
(PIT) method | Ref. DD 168], generate exponential failure 
times for the Ny copies of the first component as shown 
below, Exponential failure data is a requirement for the 
CLG and LR methods. The LG method, being non-parametric, 
can use the same data, 

a. Generate N, uniform (0,1) random numbers for 


component one. Denote these Uy31r oro ere. .» UY y 


1 
De Let Pr; 一 0 if a success occurred for compo- 
J nent one on trial j 
t if component one failed on trial 
js క N 
Appendix. E shows 
Fj = 0 if 0 <U) ; <Exp(-A;) (3,4) 
T IÍ Exp(-h;) <U,;<1 
d | OU» a o 
pt > | "P ME (ys) 
-1n(U，: ) 
e سی وت تح ات‎ 
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For the LG method, compute‏ وه 
(2N1 - 3)‏ 
1 
Ea e emi (5775‏ 
N1‏ 


and F = 2 Fy, (3.8) 


d. For the CLG method, compute the test time for 
component one, 


Tsum, = > Ts; (3.9) 
Step three: Compute 
(F) (3.10) 


for the LG method, and 


H 
A 1 
4, = -Tsum (3.11) 


for the CLG method. 

Step four: The component test times for the MLR method 
are defined differently from that of the CLG method. The 
MLR method uses only accumulated time to the last failure 
for each component; if no failures occur this is zero. 


Thus, compute for component one 


„* - ii T. irm ""o (3.12) 
1 ۱ 1 j 
j=1 
+ 
where I; is as defined on page 20. 
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step five: Compute the number of trials on component 

Two, 
SS EUREN EA 
Repeat steps two, three, and four for component two. Then 
compute Ny = Ny - [3(P, + F))] . In general, the sample size 
for the ith component is 
i-1 
E Er (5513) 

Step six: Continue alternating through steps two 

through five for the remainder of the components. 


Step seven: Obtain 


OO 7 
S= ET. (3,14) 
i=1 + 
au 
0 ers) 
m 
^ | em 
arts) = r (T;/N;) (3.16) 
jal 
A K A 
Var(Z) = E (2,/N;) (3,17) 
en 
P (3.18) 
í = 22 /var(2) (3.19) 


and the total number of replica failures 


K 
5 Mb. (3420) 


BT 
i=1 * 
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Step eight: To preclude non-integer degrees of freedom, 


denote 
[28] = max 12 smallest integer> 22) (aar) 
A A 
kan = max 12， smallest integer > 2i) (3.22) 
Step nine: For the MLR procedure, compute 
U= |min N; 1T ET =0 
" (297) 
min Tr Otherwise, 1-1,2,. Á 
and I; = rgreatest index such that 
۱ n. 
DET EU T SS (3,24) 
2. 13 1 
ae 
0 ir 8,0 
K 
amen obtain 1 =maxll, E I. (9025) 
mo 


The zero failure portion of the U computation sets U to 
the observed number of "system" successes (see Appendix D). 

Step ten: Compute the replica one LCL's for the three 
procedures, Define RD, RT, and RL for discrimination as 
below 


a. LG method 


Ra = RD = Exp جس‎ (3.26) 
9 


L(a) 


b, CLG method 


^ - -[21]2 
R = RT = Exp 一 (3327) 
Sr Dla) in a = 
c, MLR method 72 
چ‎ = RL ఇ కలా 一 (3.28) 
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Step eleven: Replicate steps two through ten 500 times 
to obtain 
RD1 ，RD2， ,RD5o0 
س9 رس‎ ۹.0 
RLy ۶٘ 020 
0 స. 


Step twelve: Order the ICL's in step eleven. 
RD(4)» స... 
RT t1)» ۳1 وه ۰۰۰۰۰۰۰۰۰۰۰ و(ج)‎ 500) 
RL ys స... 


step thirteen: Next, compute A, for the three proce- 
dures, 
RD( $00(1-a)) 
RT( $00(1-a)) 


Rb( 500(1-a)) 
and the average number of failures per replica 


500 
1 
ہد‎ కా. 


^ 
IT - 


Each procedure'!s sample variance is computed by 


500 


TO 
1 o R 
S^ EE n (AM) cd 


where Ro ula) denotes the average LCL over the 500 replicas, 
9 


Now each procedure’s accuracy can be determined by comparing 


A to true system reliability Ra» and relative dispersion 


l-a 
by comparing s, 
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Reference 1, the original LG method paper, used the 


quantity 


TT = YE N,Q. Cro) 


aS a measure of accuracy and the amount of testing necessary 
to assure sufficient accuracy of the LG procedure. As long 
as no component sample size is below ten, Table I shows 
recommended TT for the LG method (Reproduced from Ref. 1, 
at). 


TABLE I 


Recommended TT for tne Log-Gamma Method 


100(1-a)% Minimum TT Reasonable TT 
80 3 4 
90 8 10 
95 r5 20 


The simulation computes equation (3.29) for comparison 


with the above Table. 
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IV, THE INITIAL SIMULATION 


A. CASE DESCRIPTION 

The initial cases analyzed had true system reliabilities 
ranging from near .70 upwards and were composed of fifteen 
or forty components, The average number of failures per 
replica, f, ranged up to near 15. The first simulation 
used the formulas for E on page 28. The results are 
listed in Table III. Two cases were also ran with a Weibull 
failure distribution, vice exponential, to test the CLG and 
MLR procedures Sensitivity to the exponential assumption, 
Appendix E develops the PIT method for generating exponential 


and Weibull failure data. Table IV lists Weibull results, 


B. RESULTS AND ANALYSIS 

The LG method conformed to Reference 1 assertions in 
that performance improved with increasing PT, although 
Ref, 1 did not generate component sample sizes in the same 
manner as this study. This should, however, have no effect 
on the results, 

Small IT, Such as in cases one and three, resulted in 


conservative IG and CLG LOL's (ie., A, „<R.). 
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min 
max 


ہم 
LT‏ 


LG 

CLG 
MLR 
RLG 


TABLE 11 


NOMENCLATURE FOR SIMULATION RESULTS TABLES 


Sample size of 1st component; analogous to the 
number of systems tested 


Number of series connected components 
True reliability of the ith component 
True independent series system reliability 


K 
رت جا‎ 5 
5 =1 al 


Minimum number of failures in one replica 
Maximum number of failures in one replica 
Average number of failures per replica; approx- 
imation for K 

118 = N.Q; 
Sample standard deviation of Rs o) 


The 500(1-a)th order statistic of Ro La) For 
each of the three methods క 


The Log-Gamma method 
The Continuous Log-Gamma method 
The Modified Lieberman-Ross method 


The Randomized Log-Gamma method (modification 
three only) 


For example, the case one sample size for the first 


component was N, =20, and the series system modeled had 


K 215 components each with true reliability R. = „99. True 


system reliability R. = „860. There were as few as zero 


(min=0) and as many as eight (max=8) component failures 
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over the 500 replicas of case one, with the average number 
of failures per replica being 2.9 (TT = 2.9), Comparing 
A,_, to R, for accuracy, the LG 80% procedure had 099 
and s=,074, whereas the MLR hað Ay 7+851 and s 2,255, 

Thus, in case one both procedures were conservative, but 
the MLR was more accurate. However, the LG procedure had 
less than one-third the dispersion of the MLR procedure. 

Increased testing, aS in cases two and four, not only 
improved their accuracy but showed smaller sample standard 
. The same comparison holds for case five with 
case six, and cage seven with case eight, Also note the 
closeness of both Ara and s for the LG and CLG metnods. 
This basically says that even though more assumptions are 
required and information is needed for the CLG method than 
the LG method, it yields no better results, 

Case three of Ref. 1 is essentially case five of Table 
III; case five is essentially case six, and case eleven is 
essentially case two (note that Ref, 1 miscomputed TT in 
cases three and five). The accuracy achieved by case five 
was not as high as that of case three, but the sample 
Standard deviation was comparable, In the other two compar- 
able cases, both the accuracy and dispersion match closely. 

The MLR procedure is the most accurate of the three, 
both in low and high fT cases, The sample standard devia- 
tion, however, is three to five times that of the LG method, 
The most disturbing thing about the MLR method is that even 


with increased testing (Tr), s us Still high, 
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Reference 6 shows that the LR method is expected to 
have high dispersion above about ten observed failures 
(Figure 1, p. 22). The MLR method, even in the cases with 
ten or more failures, still has much higher dispersion 
than the LG and CLG procedures, Reference 6 attempted to 
improve the LR method by reduction of data loss. By 
ordering the failure times in increasing order, maximum data 
usage and chi square degrees of freedom were attained. The 
resultant ICL‘’s were found to be biased in an unpredictable 
direction (based on a Single parameter) and were even more 
widely dispersed, It was therefore proposed that all possi- 
ble orderings of the data be considered, and research is 
currently being conducted in that area. Thus this study 
pursued it no further, 

The two Weibull cases showed that the CLG and MLR methods 
were slightly sensitive to the exponential assumption, 

Since T— exponential (^) implies pi/® —. Weibull (8,+) (see 
Appendix E), then for B=,5 and T<1 mission unit, the gen- 
erated Weibull failure time is sooner than its exponential 
counterpart. Likewise, for §8=3.0, the Weibull time is 
later, In fact, case one with B=,5 gave slightly lower 
A,_, for both methods (excepting the 80% MLR ICL). With 
B=3.0, both cases gave slightly higher Aia (again, with 


the same exception as above), 
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V. MODIFICATIONS OF THE METHODS 


A. MODIFICATION ONE: EXACT DEGREES OF FREEDOM AND ADJUST- 
MENT FOR ZERO FAILURES 


1, Description 
The original LG method employed only integer de- 
grees of freedom; the initial results of Chapter IV adhered 
to this convention for all three methods. The accessibility 
of computers today makes chi square values with non-integer 
degrees of freedom readily available, Rounding up degrees 
of freedom by usage of the "smallest integer greater than 


or equal to" function can also cause monotonicity problems 


A 
159 Rs )رت‎ Consider LG 80% ICL's in the case where N,=20, 
and only one failure occurs, Then 
A K=15, 
S= سے‎ T; ae) 
pel 
= 105 
2 
r= 7% (5.2) 
å 
= 
121 و‎ 
2 
= ¿0 = 1 
165750) 
A 
then er =2 
[28]= 2 
A ^| ^ 
thus, Rs 1,(a) = Exp yi (5,3) 
1-0, 27] 





人 
Pe ae eno) | 
Es 1(.2) = LAD 12 „799 
8,2 
Now, take the exact same case and outcome, except that in- 
stead of one failure in one component, two failures occur 


in the same component. Then, 


人 


ےت 
2 
A „10 2‏ 
2.053 = = 
2f = 4,105‏ 
5 = ]= 


imeerollows that, 


R = EXD at MM 02ےے‎ 
Sec) : 7%, 
و۵‎ 
The above ICL is higher even though one more failure oc- 
curred, due to rounding up degrees of freedom 2T. 

Another improvement felt necessary was an adjustment 
for the event of no failures for the LG method, The initial 
procedure here gives 

$ c 


^ 


zal. 
and UE) OO 


Therefore, the following was incorporated: Suppose 
there are no failures in one replica with N, initial compo- 
nent trials. Then set the number of failures on the first 
component as Fj=3 and compute 2r as if one failure had oc- 
curred (this effectively constrains 2p» 2), This ad hoc and 
Simplistic adjustment will produce zero-failure LCL's as a 


function of N و‎ and satisfies monotonicity requirements. 
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The simulation was again run with the adjustments 
above for certain test cases, The results are listed in 
Table VI. ۱ 

2. Observations and Analysis 

Due to the fact that the MLR method never assumes 
non-integer degrees of freedom, it was unaffected by this 
modification. Both cases tested had slightly lower Ara 
and sample standard deviation, s for the LG and CIG methods, 


The smaller s can be attributed to less discreteness in 


2 
0 و ۰0 1 


A 
Rs L(a)* l-a 
increases faster than 2n. By rounding up the degrees of 


The lower A, , is caused by the fact that 


freedom 2n to En], tne initial simulation produced slishtly 


higher ICL's than did modification one, 


TABLE V 
Behavior of — TOR O = „2 
X l-a,2n 
ia on 
Dosen 2 
B. 





4.48 
2.97 
2,14 
1162 
1520 
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While this modification slightly decreased disper- 
Sion, it did not improve the accuracy of the LG and CLG 
methods (and had no effect on the MLR method). The only 
real gain was to solve an infrequent monotonicity problem. 
If only integer degrees of freedom are available, it is 
recommended that a more refined scheme of rounding non- 
integer outcomes be devised, For instance, in the earlier 
example, 27 = 4,105 should be rounded down to 4 vice up to 
5。 The LG ICL becomes .779, which is lower than the one- 
failure .799 outcome, but still higher than the „738 three 
failure outcome (see Table VII), 

Being that the CLG method performs no better than 
the LG method, but requires more assumptions and informa- 
tion, no further effort was expended on it. All subsequent 
modifications are directed towards improving the accuracy 


of the LG method while maintaining small dispersion. 


B. MODIFICATION TWO 
1, An Accuracy Improvement, and an Ad justment for tne 

Degrees of Freedom in the 90% Lower Confidence 

Limit Procedure 

The primary weakness of the LG method is its poor 
performance when few failures occur, More testing, and 
thereby more failures, results ina more accurate Ayia? 
Table VII lists and plots possible LG ICL's for zero through 
ten failures when N; = 20, Note the gaps caused by discrete- 
ness up to one failure. Also note that the 90% LCL for two 


and three failures is higher than that for one failure, 
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This non-monotonicity is caused by the sharp drop-off in 
سا‎ at low degrees of freedom, Ás an adjustment, 

2r for 90% LG ICL's was constrained to be at least 2.093 
the minimum 2r for 80% LCL's remains 2.0, The results can 
be seen in Table VIII. 

Consider the classical binomial component LCL pro- 
cedure on the probability of mission success p, If p is 
estimated by D=X/N where X is the observed number of 
successes in N trials, then a 100(1-a)% LCL on p, denoted 
Pra)? is the solution for p in the equation 
E (pda ا<الرو ۔‎ = (5.4 
E (;) Pp p 0 5.4) 
jx 
Under the simulation of Chapter III, if N1 is the component 
one Sample size and no failures occur ina replica, it can 
be equivalently stated that N, "systems" were tested with- 
out failure, Likewise, one component failure implies ex- 
actly one system failure, ےم ہو‎ when two or more 
different types of components fail, it cannot be ascertained 
how many system failures occurred, Figure 3 illustrates 


this for the two-failure case, 


Figure 3 
Component Failures Versus System Failures 
Component 1 2 K 1 2 X 
N, systems ONSE wee = N, systems Be. > 
tésted with _ _ _ op tésted with _ , l 
one system K En ۴ two system : ELE 
failure 5 failures ? : 


5 


Treating a "system" as a component, equation (5.4) 

can be used to calculate Rg een) when zero or one failures 
9 

occur. 

a. Zero failuress 

For X=N successes, D 1ی ے‎ Modifica- 
: PL, (a) ; 

tion two revised this upwards to 


Rs, 1(0) = + 21 - ja) (555), 


b, One failure; 


Equation (5.4) was applied directly for X=N,-1. 


A 
That is, R is the solution for p in the equation 
S,L(«) 


N , : 
2 (Mo) pd(t-p)Nt-d = q (5.6) 
J=N, -1 J 
or N, p 171 (1p) * pl = a (5:7) 


C. Two through eight failures: 


The LG method was adjusted upwards as follows: 


Rs L(a) E رو ری وخ‎ 3 Rs 1(0) 0 + PTA | (5.8) 


A 
where R (N,,FT) denotes the LG LCL computed with the 
۳/0 ( ۰ 1 
observed total replica failures FT when the original sample 
size ís Ny » 
d. More than eight failuresi 
The original LG procedure is used with exact 


degrees of freedom, The results are in Table IX. 
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2, Observations and Analysis 


Modification two did improve the accuracy of the 
LG method, while slightly increasing dispersion. The 90% 
procedure was noticeably improved by the degrees of free- 
dom constraint, although still conservative, There is, 
however, Still too much discreteness in values of i 
when few failures occur, Table X generated at NPS by APL 
Subroutine "Histlist" (see documentation), illustrates this 
for case five, This discreteness usually causes the proce- 
dure to yield more than the target l-a confidence, and is 
a problem with nearly any method using only attributes data. 
This being the case, the next logical step was to 


fit a continuous curve to the procedure by randomization, 


C. MODIFICATION THREE: RANDOMIZATION 
1, Randomization of the Log-Gamma Method 

Many LOL procedures which use only attributes data 
are degraded in that they yield conservative bounds, Be- 
cause only a finite set of outcomes can be observed for any 
test plan, then only a finite set of LCL's can be realized, 
Thus, e is a discrete random variable which ideally 
gives the ICL closest to but not greater than the "exact" 
LCL. 

The 1975 article of Barr and Jayachandran | Ref. و7‎ 
D 67-68] developed an exact, randomized binomial component 
LCL procedure. The same rationale used before (see Figure 3) 


allows the randomized binomial procedure to be used for 


48 





TABLE X 


Case 5 Log-Gamma 80 Percent ICL's with Modification Two 


SERIAL NUMBER 


ORDERED DATA 


„407610 
اباب‎ 90 
„tl:960 
"71830 
. 475840 
177390 
582700 
. 484660 
. 487430 
. 527610 
. 534450 
. 538120 
: 92 ٤ 
. 540650 
. 541090 
. 542420 
. 542860 
.571020 
. 576210 
TAS 
576500 
. 580530 
, 581070 
5-5 950 
. 584260 
. 584910 
¿616160 
«620600 
» 622570 
60ء‎ 
. 624390 
„625090 
„627950 
„628010 
.628080 
. 662770 
. 666380 
. 666490 
„667840 
„668570 
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OCCURENCES 


p 


pa 


¡a 


3* 3€ 


pa 


RH 


pa 


HHX 


N 
OPEFEPRIWLHDURRKRRNFÆFPRFVYOYØØANVNRHONVNVPRRRDNPRPPRPROVWUPRPARRNUNVDPRAP RR 


* 3+ 3¢ 


۲ 


*% 


pa 


چد +3 


pa 


PERCENT 


„002 
„002 
„002 
„002 
„004 
„ 006 
„004 
„002 
„002 
„002 
„ O06 
„020 
„002 
o 002 
„002 
‚004 
„002 
„022 
„002 
„GIL 
. 00۱ 
„020 
„002 
, 002 
‚00% 
. 004 
„026 
„010 
„028 
„002 
„ 048 
„ 004 
„002 
„002 
„002 
, 05 
„002 
„028 
„002 
„032 





SERIAL NUMBER 


188 
209 
210 
216 
218 
219 
220 
257 
279 
304 
310 
312 
0ئ‎ 
380 
01 
403 
m 
469 
490 


TABLE X 
(continued) 


„669380 
„669530 
„670340 
„671090 
„672070 
. 673220 
. 710230 
2250020 
, 714870 
„715380 
„716280 
. 756600 
755650 
.759200 
-761320 
. 789380 
. 791340 
. 890140 
‚938140 


ORDERED DATA 


50 


OCCURENCES 


21 


3 3t 3* 


* 


e HH 
3 + + 

* 3t 
* 


H H X H H Ý 
+ A ¢ 
3e 3€ + 


3t 3€ 3€ 3€ 9€ 
3e 3t 3t 

A H 4 

* 


PERCENT 


042 „ 
290027 
012 
004 
002„ 
002 
با07 , 
بابام , 
050„ 
012„ 
004„ 
باون , 
e 042‏ 
042 , 
004 , 
078 
054 ,+ 
042 , 
ULE‏ 





"system" ICL's when zero or one failure occurs, Equations 
for the Randomized Log-Gamma method (RLG) are derived in 
Appendix F. An example is given in Appendix G, 
a. Randomized Log-Gamma Procedure 
(1) If zero failures occur, generate a uniform 


(0,1) random number Y, Then 


1/N1 


R లంట 
Rs,L(a) = (Im if Y<l-a 
(5.9) 
1,0 if Y>1-a 


(2) When one failure is observed, generate a 
uniform (0,1) random number Y, Then Ra D is the solu- 
9 
tion for p in the equation 


N1-1 N1 


(1-Y)Np (1L-p)+ p * ہے‎ if Y<l-a (5,10) 


This solution can be approximated by 


Ra نا‎ ۰۳ 


for Y>1-a (see Appendix F). 
(3) Two through four failures inclusive: 
Generate a uniform (0,1) random number Y; 
for each component that nad at least one failure, i-1,2,... 
...»K. Then use the LG method with P =F, -Y instead of 
Fs. Thus, the RLG formulas become 
2N, - 3 


Ay = و211‎ - 2 (5,12) 


B; = gg (RS) 
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— i i ET i i 
ہت و تس‎ ME) (5,14) 
i i 
A 97 
S = 7 Ti 65815) 
1= 
^2 
E ہے‎ (5.16) 
Se 
qe Ue 
and E oid xD — నా (5.17) 
1-a,2r 


The only differences between the original LG formulas and 
those above is the computation of i in equation (5.154), 
and usage of exact degrees of freedom, 
(4) Five or more failures: 
Use the non-randomized LG method with ex- 
act degrees of freedom. That is, use the formulas in (3) 
above except use Fi vice ED in the computation of ie 
2. Modification Three Cases 
The number of cases Simulated under modification 
three was expanded greatly. Cases one through twelve have 
an N, of either 20 or 50, and cover the range of pertinent 
system reliability Rye All twelve cases were first ran 
with the number of components K=15, and then reran with 
K=40, Cases with equal component reliabilities were also 
ran with K=4, The system sample size was reduced to ten, 
and the applicable twelve cases ran again, with K=15, 40, 


and 4, Modification three results are listed in Table XI 
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య క గక ౭౦. ఆసన for K=40, and Table XIII for K=4, 
Table XIV lists results of cases from Ref, 1, which had 
highly-varied component sample sizes, 

3. Results, Analysis, and Conclusions 

All three metnod results were tabled for comparison 
under modification three, although randomization of the LG 
method was the only change. 

The accuracy of the RLG was comparable to that of 
the MLR, while having much lower sample standard deviations, 
Both the RLG and CLG methods had dispersion inversely pro- 
portional to sample size, while the MLR continued to con- 
sistantly have high dispersion regardless of the circumstances, 

The most significant improvement in the RLG was its 
performance with few failures and less testing, Especially 
in the K=4 and K=15 cases with small T, the RLG was near 
exact while the CLG sometimes gave A, , of 1,00 (meaning 
at least the upper (100)a% of the simulated distribution 
was ICL's of 1,00). The MLR procedure was inaccurate for 
very high reliability systems, such as case 11 with K=15, 
due to the upper limit of its possible ICL values, In 
case 11, with N, =20, the maximum MLR LOL is .923 with zero 
failures (see Appendix D), whereas true system reliability 
is „956. As in the earlier LG results, the 90% RLG procedure 
was usually less accurate, less stable and slightly more 


dispersed. The presence of one component with dissimilar 


8 





reliability, like case seven with K=15 seemed to degrade 
accuracy somewhat, Additional testing, however, as in case 
eisht, dampens this out, 

Cases one and two of the original LG paper [Ref. 1] 
were also simulated with modification three to check sensi- 
tivity to mixed sample sizes (see Table XIV). Note that 
randomization was rarely necessary in case two (min=2), 
and never in case one (min=5). The accuracy of all three 
methods was high due to the large amount of testing. Note 
the small sample size of the twelth component in case one, 
This seemingly had little effect on the RLG method's accura- 
cy, although s is suspiciously large when compared to case 


two, When R was set to .85, the effect on the RLG became 


12 
clearer, When a component twelve failure occurred, its 
contribution to S and Ë so dominates that Rói cu) is nearly 
halved from what it would have been otherwise. 

Figures 4 and 5 give histograms for Rs (a) of the 
RLG procedure generated under cases one and two respectively. 
Both have skewness and kurtosis roughly comparable to a nor- 
mal distribution, and both have medians near-equal to their 
medians (symmetric). Figure 5 (case two with N, = 50) has 
a lower maximum 0 and a higher minimum than Figure ہا‎ 


(case one with N, = 20), reflecting the decreased dispersion 


of the RLG procedure with increased testing. 
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APPENDIX A 


APPROXIMATION OF S BY A TAYLOR SERIES 
The Taylor series approximation of f(x) about x=a 


can be stated as: 


R 2 
f(x) -f(a)* f'(a)(x-a) ايقل‎ + 
یب عتتقتے‎ (4, 1) 


Here, f(x)=-ln(x), where x=(1-Q;¡). Let a=1， 
Then (x-a)= -Qi 
Ma) =0 
۱ T: uw. ' ES 
I (x) 1-0 » f (a) = l 


fer) سس‎ Sa! 
(1 -Q;) 
f'(x) = ——À^ f"(a) = 2 


Equation (A.1) becomes: 


ر 2 
t‏ اع ,انالا ,10( +0= f(x)‏ 


2 . 
í 5 E 
Biome) We ete ۴۳6 دب‎ 
ص‎ 3 
= E (22) 
jz1 ( 
K ×۰ کا‎ Q; 
Thus, 5 = JE -In(1-Q,) کے‎ E 2 GE (A.3) 
i=1 i=1 j=1 2 
Dropping terms of order three and higher, 
: A (4.4) 
S یم‎ às Q. + چ‎ A, 
¡=1 | it 2 | 
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APPENDIX B 


AN UNBIASED ESTIMATOR FOR Ti 
2 
人 
Let Ti = Qi +, An unbiased estimator I. is desired 


for Ti, which is of the form 
B 


T. - AQQ iQ? (B.1 
; = Aiai + 25 Dei) 
A PF. 
where ai = 天 , The number of failures of ith component F; 
i 
is a binomial (N = Nes p=Q,) random variable with mean N.Q. 
and variance N,Q; (1-9; yar Thus. 
q — E(F,) = (3.2) 
E(Q1) = N; (F; ) = qd; BZ 
Q.(1-Q.) 
5 =- 二 
Var(Q,) = 2 Var(F,) = గా (B.3) 
i 
i 
A A A 
E(QÍ) = Var(Q,) + (E(Q,))* (3.4) 
= ELT 5 + a“ 
Ni 1 
2 2 
۰۰٦ 
s 
Substituting 


^ A A 


Bs B, 5 
1 1 
A 
That Ti is unbiased implies 
Qe 
E ఆఆ - adi. 
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Equating coefficients with equation (B.5) 


Bi 
1 
B. 
జు _ 
B. = N; = 1 (B.8) 


Solving the above two equations simultaneously gives the 


desired expressions for Ai and B; 


EN; - 3 
As E (3.9) 
1 
Ð 
B; سوہ وو‎ (3.10) 


2 
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APPENDIX C 


DERIVATION OF THE ACTUAL VARIANCE OF E 


K T: 
AND COMPARISON WITH C = 
i-1 N; 
Reference 1 used T as a variance estimate for Es 
1 
Then for the method of moments gamma fit, 
A 
A K 7 TE 
Var(S) = 2 Var(T, ) = 2 X 
i=l jT 


|- 


(CL) 


The following is a derivation of the actual variance 


^ 
of T., and then a comparison with the above estimator. 


y 
B. 
: NN 2 A 
Since T: - AQ, + 3 Q; : (62+ 
B. 
^ » 2 ۸ A رھ‎ 7 ^2 
then Var(Ti) - A;Var(Qi) * A,B,COV(Q. ,Q7) * T Yar(Q?) (1635) 
A Fi 
where کو ہوک‎ The total number of failures on the ith 
d 


component Pi, is a binomial (N=N,, p=Q,) random variable, 


Then 


E(Q,) = E(F,)=Q, (C.4) 
l 
Q,(1-q,) 
A 1 1 
Var(Q, ) = vé Var(F. ) = a TE (655) 
E(@2) = var(¢ a, ))“ 
2) = Var(Q,) + (E(Q,)) 
Q, - QF 
i 2 
= q; (C6) 
L 
and cov(ä,,02) = E(Q3) - B(Q, )E(Q2) (0.7) 
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To solve equation (C.7), E(Q7) 三 Lu must be solved, 
INS 
i 


Model I = | 1 if the mîn test fails 
0 if the men test does not fail 
N; 

and let Fi رح‎ In (Gro) 
m=1 

Then F? = ےج ےہ‎ 

and E(FS) = sya) E(Ij) + rr E(15)E(1,,) + 


(all three equal) (two the same and 
one distinct) 


۳ - 
0 ۶ ۹۶۰۶ 2س‎ (C.10) 
(all three distinct) 


2 M 2 3 
= N,Q,+ 3N, (N, - 190; + N.(N, - 1)(N; - 2)Q1 
This implies E(Q?) = 3 E(72) 
_ k E 2 _ ø 5 
- ۲ + 3(N; - 1)QÍ + (N; - 1)(N; að 


T 
Now, using equation (C.7), 


N; 
1 


A Q. 
cov (Â d£) = 号: + any 38, 200, nt | VE 


All that remains unknown in equation (C.3) is Var(Q“). 


E(Q;) - (E(Q2))? 


1 4 Qi - 7 2 Á (C.12) 
HUE ).- DX OR. 
N; i 


4l 


A 
Var(QÍ) 








Using the same model as before, 


z E(IÍ) Ber 


x E(12) E(14)‏ تر مو 


E(F; ) = 2)5 ) + 


5 k 


(all indexes ns equal bes indexes the 
equal) in pairs) same and one distinct) 


SE E 2 
۲1 E(T,) E(1,) E(14) 


(two indexes the same and two distinct) 


E 


ESLE 
Iek mn 


+ 


E(1,) E(I4) E(1,) E(Ipp) (0.13) 


(all four indexes distinct) 


2 2 3 
N,Q, + IN; (N; - 1)Qj + YN; (N; - 1)Q; + 6N; (N; - 1)(N, - 295 
۷ 
+ N.(N; - 1)(N, - 2)(N; - 3)Q, 
Exon (0.13) is used to find 


- 
(గి) = " + 7(N, - 1Q + 6(N, - 1) (N, - 2)03 


+ (n, 1008, - 2900, = 3198 | (C,14) 


Substituting equation (C.14) into equation (C.12), 


zZ [| శు 
= Vl: 


: zi = DN = 2962 | (C.15) 


The final result can now be obtained by substituting equa- 
tions 615), (C.11), and (C.5) into eduation (0.3), yielaihs 


TE 





۰ 0 


Var(T.) کت‎ Eus 0 14 (2N. - 3)Q. - 2(N, - 1)Q* 
i N; i i a ine 1 1 1 1 
1 
Bí Qi 2 
Eus 3 1+ (ÓN, - 7)Q; + (N; > 1) (N; - 3) ( 
1 
2 
- 2(N. -1)(2N. - | (C.16) 
1. 1 i 


The cumbersome expression above can be greatly simplified 


by dropping terms with divisors of NE and higher, yielding 


Var(T,) = Maeg (A; + B. TH ) d (C717) 
N; E 
SS ME ^ 2 Sg 
Then Var(T. ) (1-4. ) (A; +B. ¿9 ) | where Q. im (6515) 
N; L 


A comparison was made between Tar(î, ) as given by equation 

(C.18), and that used by Ref. 1. The results are listed 

por N4 =O, 2O ana 50 in Table XV. The simplistic esti- 
K T, 


mator 2 M turns out to be an excellent approximation 


isl í 


^ K 
for Var(S) - A Var(T, Ye 
=1 


aÐ 


TABLE XV 
Actual Variance of T. Compared to Estimate 
(One Failure Per Component) 
K T. ۵ 


A N JA KQ 
FAILURES var(8)= y = Var(S)= E TO -â )(A,+B, 7 ye 
اد‎ 


وآ 1ے و 

1 N,=10 „01000 „01003 
20 „00250 „00250 

50 „00040 . 00040 

2 N,=10 5 ‚02242 
20 „00527 „00527 

50 . 00082 „ 00082 

3 N,=10 „03469 . 03480 
20 . 00804 .00805 

50 90123 . 00123 

4 N,=10 . 05032 . 05050 
20 „01113 „01113 

50 . 00167 „ 00167 

5 N,=10 . 06594 . 06619 
20 ‚01421 ‚01422 

50 ,00210 ,00210 

6 N,=10 ,08635 . 08672 
20 01767 01769 

50 „00255 „00255 

7 N,=10 . 10676 „10725 
20 ro TIS . 02115 

50 „00301 „00301 

8 N,=10 213454 9526 
| 20 „02504 „02506 
50 . 00348 , 00348 

9 N,=10 ‚16231 „16327 
20 „02895 „02897 

50 „00395 „00395 

10 N,=10 „20231 ,20377 
20 „03339 „03342 

50 „ 00445 , 00545 
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APPENDIX D 
DERIVATION OF THE ZERO FAILURE CASE 
FOR THE MODIFIED LIEBERMAN-ROSS METHOD 
The variables U and I are undefined when no failures 
occur, Consider type III censoring in which N components 
are placed on Simultaneous test, with replacement, until a 
pre-selected time To [Ref. D, 238]. Type III censoring 
is a component ICL procedure, and not a system ICL procedure, 
Now, consider the Modified Lieverman-Ross (MLR) procedure 
where Ni items of the ith component have been tested with no 
nes, 1=1,2,,s00K. If the N; are identical (as is the 
case in this study for zero failures), then the above is 
equivalent to testing N; "syStems" with no failures until 
م1‎ < 1 "system" mission. 
The type III 100(1-a)% ICL on the mission probability of 


Success p of a component is 


2 
Pra) = EXP (sena (D.1) 


where F is the observed number of failures, N is the number 
of components on simultaneous test, and To is the preselect- 
ed test time, Then for F=0, 2F+2=2 and 

2NT, - 2N, " 1 system mission unit - 2N; 


A 100(1-a)% ICL on system reliability R¿ When no failures 


occur is ^ E Ai (D,2) 
sia), = SP EN; 


for Ni edual, ول‎ KR。 If the Ni are not equal, the 


minimum Ni should be used. 


75 





APPENDIX E 
USE OF THE PROBABILITY INTEGRAL TRANSFORM (PIT) IN 
GENERATING EXPONENTIAL AND WEIBULL FAILURE DATA 
Assume an exponential failure process for the ith com- 
ponent with rate Mi per mission unit, Then the reliability 
of the ith component for completing one mission unit isi: 
一 E =æ = a 9 
R¿(1) xp ( ^s 1) (F,1) 
Utilizing the PIT method, draw a uniform (0,1) random num- 


ber, denoted "U", and form the equality: 


U = Exp (-\, +) (2 
Solving for time-to-failure t 
ln(U) s “r,t (F,3) 
It follows thats 
t = =U (Fb) 
D 


The computer program uses the following algorithm: 

1) Generate a uniform (0,1) random variable 

2) is U = Exp( -*.t) is less than Exp(-h, ° 1), then t>1 
mission unit, and a success has occurred. Truncate the 
time-on-test for that component trial to 1, 

3) Otherwise, a failure has occurred at t = in) 

The simulation was also ran under a Weibull తతా! 


process for sensitivity analysis. The Weibull (B,)) densi- 


క్‌ m4 
ty function is بو = رت‎ Exp( -A£P ) (F.5) 
and hazard function 

n(t) = prt? t>0,%>0,8>0 గర) 
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Using the PIT method as before, draw a uniform (0,1) 


random number U and set 


U = F(t) 

= Exp(-h, +") (E7) 

thus T = (ip) ‚Up (F.8) 
al 


Comparing equations (F.4) and (F.8) shows that if 


y == ni is an exponential (A; ) component failure time, 
1 


then pi /P is Weibull (8, Mi). Thus, to perform the Weibull 
sensitivity analysis, only the above change needs to be 
made to the simulation algorithm (Chapter III). 

Taking the derivative of the Weibull hazard function 


with resvect to time 


AE = స (0 قرو‎ (F.9) 


dt 
Therefore, any O-«8-«1 yields a decreasing hazard rate, and 
any B>1 gives an increasing one, Case one was ran with 
B=.5 and 3.0. Note that the comparison logic, used by 
the simulation algorithm to determine if a failure has 
occurred, reduces to the original expression. The number 
of failures is then the same as before. The Weibull results 


are shown in Table IV. 
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APPENDIX F 


DERIVATION OF THE RANDOMIZED LOG-GAMMA METHOD 

The technique for randomization of the classical bi- 
nomial model is a direct application of Barr and Jayachand- 
ran's 1975 article | ref. 7|. It can be applied directly 
to obtain a system reliability ICL when a total of zero 
or one failure occurs among all component tests, because 
tnese cases translate into equivalent system failures (see 
Figure 3, p. 45). 

The binomial point estimate for p is p-X/N, where X is 
the number of observed successes in N trials. The classi- 
cal binomial 100(1-q)% ICL for p, denoted P (a) is the 
Solution for p in the equation 
۶ دی رای‎ ) da (F.1) 
j=x * 
when no failures occur, or X=N, this becomes 


a‏ ا ره 


and for one failure, Pr (d) is the solution for p in 


Np (1 -p)+p =a (F,3) 


This Dr (a) has the property that Pr MET >1-0 
Reference 7 develops a method of randomizing discrete 

LOL procedures which are conservative, and illustrates it 

for the binomial case, The exact randomized 100(1 -a)% 


binomial LCL for p is developed as follows. 
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a. The number of observed successes X is a binomial 
random variable with probability mass function B(x;N,p). 
Define h(p) as the smallest integer such that ° 

Pr(X> h(p)) <a o, ا سس‎ (F.4) 


Thus, h(p) denotes an observed value of X. 


b, Define 
a= Pr(X>h(p)) (F,5) 
and roa rt = (a= a BMP) IND) (F,6) 


where F is the "randomizing" distribution function of Y. 
Reference 7 shows that F can be the uniform (0,1) distribu- 
tion function without loss of generality. Then using 
equation (F.5), y(p) can be simplified to 

y(p) =1-(a-Pr(X>h(p))/Pr(X=h(p)) (F.7) 
Finally, define 

g(p) = h(p) + y(p) (E78) 

Then the exact binomial randomized LCL for p, L (x+y) 


is given by 
+ g (x y); 1 - a < x+y <S N+1-«a 
L (x+ y) = 
0 ; x*y «1-a (F.9) 
1,0 ; x+y > N+1-0 
To solve for L’(x+y), forn the equality 


g(p) = X+Y = h(p) + y(p) (9,10) 
Since h(p) is an observed value of X, this reduces to 


Y = y(p) =1 - (a = Pr(X>h(p)))/Pr(X =h(p)) 65011 ( 


(2 





By solving for p and applying equation (F.9) 


the solution for p in 


+ 
L (x + y) = 
(1 - Y)Pr(X=x) + Pr(X>x)=a 1-a < x+y <N+1-« 
(F,12) 
1,0 x+y > N+1-a 


+ 
Thus, for no failures, or X=N, L (x+y) is the solution 





OP pin 
(1-Y)p" =a (F,13) 
R a y 1/N y<t-a 
in 0+۷ 8 ۹۹۷, 5 Bon 
170 y>1-c0 ` 


* 
when X=N-i, or one failure, L (x+y) is the solution for 
pin 
MEY) EAX =x) క x)= a 


(F.15) 
Or లం అ ల « 


When two or more failures occur, a component ICL pro- 
cedure such as L'(x-* y) cannot be used, The original Log- 
Gamma procedure can, however, be "randomized" as follows. 
Draw a uniform (0,1) random number for each type of compo- 
nent that fails. Denote it Y; for the ith component, and 
compute the "randomized” number of failures 


* 


F = F; - 1; (iði, 16 ) 


where F; is the observed number of failures on the ith 


component, 
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The Randomized Log-Gamma procedure (RLG) was developed 
by trial and error via computer simulation, using the above 
adjustments where most effective, It can be stated as 
follows: 

Randomized Log-Gamma Method 


a. Zero failures in all components: 





2 1/N 
R Ü_ | (F=) ysı-a 
1.0 ۳ < 1 ه‎ ర 
where Y<a uniform (0,1). 
b. One failure among all components: 

7 the solution for p in 

R = (F.18) 
Sala) N1 -1 Ny 

(1-Y)Np (nb ries = a 


The use of 1ا‎ as an approximation to వ oy ftor. y < و‎ 
y 
1s reasonably accurate for low N, and very accurate for 


large Na» as shown below, 


TABLE XVI 
Accuracy of the Approximation oM 
و ےت‎ 
2 
A 
Actual Sad) No = 10 20 50 
UN „837 „915 „965 
051 2925 „968 


Co Two through four component failures: 
For each different type of component that had at 


least one failure, draw a uniform (0,1) random number, 
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Denote it Y; for the ith component, and compute 


F. =, - Y. (F.19) 


Compute Log-Gamma method with exact degrees of freedom and 
P; instead of the observed number of failures Pi. 
d. Five or more failures: 
Use the Log-Gamma method with exact degrees of 


freedom, 
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APPENDIX G 


AN EXAMPLE OF THE RANDOMIZED LOG-GAMMA PROCEDURE 
Let N, =20, K=15, and a=.2. Suppose that the replica 
has no failures, Then a uniform (0,1) random number Y is 


generated; say Y-2,5. Then since Y«1-0, 
R 2 
s .L(.2) 7 (TF) = +755 


as opposed to the classical binomial ICL of .923. 

Suppose some future replica had one failure in some 
component, and Y=.6 is RM Then Rs 1.2) is the 
solution for p in the equation 


20p:?(1-5)(1-.6) 4 p*? = .2 


Therefore, ETE = „092, The classical binomial ICL is 
„858. 
Now, suppose another replica had two failures in the 
5th component and one failure in the 15th component, and 
i 


that و‎ 2 and f15= 8 were generated, Thus, F„=1.7 and 


* e 
F15 =，2。 Bp 0 LG procedure, N¿= 20, Ns 719 


Ss = 0 9 一 „0969 
^. 2 
۲ SEEN MOSS i 
Be As ER 6 
تب رد‎ 
ii 
2r = 3.853 


2 ۲ 
[26,3853 * 1:56 
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Then E Be -22S 
Rs L( .2) = Exp m 
2 


without randomization, the ICL would have been 


Rs 4 
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APPENDIX H 


THE COMPUTER PROGRAM 

The following pages list the computer program used for 
the simulations. Flags were used to designate the portions 
of the program to be executed, depending on which modifica- 
tions were desired, Proprietary subroutines [Ref, 9 | were 
utilized for inverse chi square values (MDCHI), sorting 
the approximate distributions of Rs aa)? (VSORTA), gener- 
ation of uniform (0,1) variates (GGUB), and solving for 
the RLG procedure (ZREAL2). The APL subroutines HISTLIST 
and HIST were used in Table X and Figures 4 and 5 respec- 
tively [Ref, 10 | . The plotting subroutine PLOTP was used 
in Tables VII and VIII [Ref. 11]. 

The program was written in the FORTRAN IV language for 
use on the IBM-360 system at the Naval Postgraduate School, 
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